All-optical memory based on the injection locking bistability of a two-color laser diode 
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We study the injection locking bistability of a specially engineered two-color semiconductor Fabry- 
Perot laser. Oscillation in the uninjected primary mode leads to a bistability of single mode and 
two-color equilibria. With pulsed modulation of the injected power we demonstrate an all-optical 
memory element based on this bistability, where the uninjected primary mode is switched with 
35 dB intensity contrast. Using experimental and theoretical analysis, we describe the associated 
bifurcation structure, which is not found in single mode systems with optical injection. 

PACS numbers: 42.55.Px, 42.65.Sf, 85.60.Bt 



INTRODUCTION 

Much current research in optical device physics is fo- 
cused on bistable, switchable semiconductor lasing ele- 
ments. The aim is to provide sophisticated functions 
for all-optical signal processing, including optical mem- 
ory, flip-flop or optical logic [1-4]. The systems consid- 
ered comprise single [3| or coupled [I!, 0] devices includ- 
ing micro-cavity lasers Q, edge-emitting devices [t^ and 
VCSELs [8-10]. Although laser diodes are known to ex- 
hibit complex dynamical features, they are nevertheless 
successfully modeled by surprisingly low-dimensional and 
generic systems of equations. Laser diodes therefore pro- 
vide model systems for the understanding of fundamen- 
tal concepts of nonlinear dynamics, such as non- invasive 
chaos control [ll| or synchronisation of coupled chaotic 
systems [l^ . 

Switching between co-existing states of bistable semi- 
conductor laser devices can be achieved using modulated 
or pulsed optical injection. The structure of the bifur- 
cations in single mode semiconductor lasers with opti- 
cal injection are well understood. For a recent review 
see reference [l3]- While two-mode models have success- 
fully explained the switching between single mode states 
with different polarizations in VCSELs 1^ , a general un- 
derstanding of nonlinear dynamical phenomena in laser 
diodes for the case of more than one simultaneously ac- 
tive mode is still lacking. In edge-emitting lasers a rea- 
son for this is that plain edge-emitting Fabry-Perot (FP) 
laser diodes will generally operate with many longitu- 
dinal modes close to threshold, which complicates the 
interpretation and modeling of experimental results con- 
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Using an inverse problem solution 
the threshold gain modulation in wavenumber space to 
the spatially varying refractive index in a FP laser, we 
have designed a two-color FP laser diode with THz pri- 
mary mode spacing. The large engineered mode spacing 
of this two-color device leads to weakly coupled modes 
with the result that the bias current can be adjusted 
so that the two primary modes oscillate simultaneously 
with the same average power level. With only two active 



modes, such a device represents a test-bed for the system- 
atic study of the nonlinear dynamics of multiwavelength 
semiconductor lasers. 

In this paper, we examine the locking bistability of the 
two-color device with optical injection. In contrast to the 
case of a single mode laser, where a bistability is found 
between a stable limit cycle associated with wave mixing 
and the injection locked state 13, IQ], we demonstrate 
that in the two-color laser a bistability between the in- 
jection locked (single mode) state and a two-color steady 
state is possible. We experimentally demonstrate all- 
optical switching between these bistable states, where the 
uninjected mode of the device performs an inverted mem- 
ory function with greater than 35 dB intensity contrast. 
We explain the underlying mechanism using modeling 
and bifurcation analysis, showing that a four-dimensional 
model of the system provides excellent agreement with 
experimental results. Our results demonstrate a new and 
unique approach to in-plane all-optical signal processing. 



EXPERIMENTAL RESULTS 

The device we consider is a multi-quantum well 
InP/Ina;Gaj,Ali_:c_yAs FP laser of length 350 /im with a 
peak emission near 1.3/im. It incorporates slotted regions 
etched into the laser ridge waveguide that determine the 
lasing mode spectrum. Further details of the design and 
free running lasing characteristics of the device can be 
found in [l9|. Adjusting the device current, we equalize 
the time averaged optical power in each primary mode 
of the free running laser. The primary mode spacing be- 
tween the short wavelength mode vi and the long wave- 
length mode V2 is chosen as four fundamental FP modes 
(480 GHz). 

For our experiment, we use the two-colour laser in a 
master-slave configuration, where the master laser is a 
tunable laser with < 100 kHz linewidth. A hysteresis 
loop of the light output power in mode vi for increasing 
and decreasing injected power in mode V2 is shown in 
Fig. [H In the insets of Fig. [1] we show the co-existing 
states associated with the locking bistability. The in- 
jected mode V2 is at long wavelength and the detuning 
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FIG. 1: The optical power of the short wavelength mode, v\, 
as a function of the injected power at the long wavelength 
mode, i^2, in up and down sweep as indicated by the arrows. 
The lower left inset shows the two color equilibrium state. 
The upper left inset shows the injected mode in this state in 
close-up where a two-peak structure is seen. These peaks are 
separated by the detuning frequency. The right inset shows 
the single mode locked state of the device. 



in this case is Aw = — 14 GHz, wliich is approximately 
3.5 times the relaxation oscillation frequency, i^Ro = 3.9 
GHz. For weak injection the dynamics are two mode 
and time-periodic due to the wave-mixing response. At 
large injection, single mode oscillation with locking of 
the output to the injected light is observed. Decreas- 
ing the optical injection from inside the locking region, 
the single mode locking state and a two-color state are 
found to co-exist. In the two-color state we find that 
the uninjected mode has a greater intensity than the in- 
jected mode, whose spectrum (upper left inset of Fig. 
[ij comprises two distinct peaks, which are separated by 
the detuning of the injected light frequency from the free 
running cavity mode frequency. 

The optical spectra in Fig. [1] are complemented by 
the mode resolved power spectral densities depicted in 
Fig. [2l We observe four regions of qualitatively different 
dynamics, the boundaries between which we have indi- 
cated using corresponding values of injected field strength 
Ki, i = 1,2,3. At weak injected power (injected field 
strength K < Ki) we find a two- mode wave- mixing re- 
gion where the detuning frequency is observed in the dy- 
namics of the injected mode and a weak feature at the re- 
laxation oscillation frequency is observed in the dynamics 
of the uninjected mode. With increasing injected power 
{Ki < K < Ks) this feature in the uninjected mode 
becomes considerably stronger and at the same time the 
peak at the detuning frequency in the dynamics of the in- 
jected mode broadens. At large injected power {K > K3) 
we enter the single mode locking region where the spec- 
trum is structureless. For K2 < K < K3 we identify a 
region of optical bistability between a two-mode and a 
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FIG. 2: Power spectral densities (PSD) of each of the pri- 
mary modes of the device as a function of the injected power. 
The upper panel shows the PSD for increasing injection, while 
the lower panels show the PSD for decreasing injection. The 
frequency detuning is Auj = -14 GHz and the relaxation os- 
cillation frequency is i/^o — 3.9 GHz. Left panels: uninjected 
mode, ui. Right panels; Injected mode, 1/2- 



single mode state. The selected stable state depends on 
the sweeping direction. 

We have been able to switch between these two bistable 
states using an appropriate pulsed modulation of the in- 
jected signal. The experimental setup and the mode re- 
solved time traces in response to optical injection mod- 
ulation with 5 ns pulse width are shown in Fig. [3l One 
can see that the device operates as a memory element 
where we can regard the mode at short wavelength as per- 
forming the (inverted) memory function and the injected 
mode at long wavelength as the writing channel. In par- 
ticular, the intensity of the uninjected mode is switched 
between an "off" state and an "on" lasing state with a 
contrast ratio of greater than 35 dB. 



MODELING AND BIFURCATION ANALYSIS 

To understand the structure of the bifurcations that 
lead to this bistability we have adapted the well known 
model of a single mode laser j^^l with optical injection 
to account for the presence of a second lasing mode. The 
system of equations in normalized units may be written 
as follows: 

Ei^ ^{l+ia)igi{2n+l)-l)Ei (1) 
E2 = [^(1 + «a)(52(2n + 1) - 1) - iAuj] E2+K (2) 
Th = P-n-{l + 2n)Y.^^gm\Era? (3) 

where the nonlinear modal gain is 



l + e^/3,„„|£;„| 



(4) 
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FIG. 3: Left panel: The experimental setup schematic for the 
all-optical memory element device operation; PC, polarization 
controller; A.M., amplitude modulator; BS, beam splitter; L, 
lens; I, isolator; G, grating. Right panels: The intensity time 
traces of the uninjected mode, ui, the injected mode, 1^2, and 
of the injected field strength. The pulse duration is 5 ns. 



In the above equations, g^m is the linear modal gain 
while the ePmn determine the cross and self saturation. 
The normalized complex electric field amplitudes of the 
modes are given by Em and the normalized excess carrier 
density is n. T = "/Ts, where Ts is the carrier lifetime and 
7 is the cavity decay rate. The phase-amplitude coupling 
is given by a while P is the normalized pump current. 
The bifurcation parameters are the normalized injected 
field strength K and the angular frequency detuning Aw. 

Parameters used are a — 2.6, P — 0.25, (50 % above 
threshold), = 0.00125, e = 0.01, /3i2 = P21 = 2/3 
and (3ii = (322 — 1- This ratio between the cross and 
self saturation is consistent with the stability of the two- 
mode solution in the free running laser [l9| . 

Note that although we have provided a complex equa- 
tion for the field Ei , the phase of Ei is in fact decoupled 
leading to a four-dimensional system of equations. Only 
the intensity of the uninjected mode infiuences the dy- 
namics, reflecting the fact that the mode spacing is in the 
highly non-degenerate regime. More generally, our model 
describes two oscillators, whose amplitudes are coupled 
via a global variable n, and where one of the oscillators 
is forced by a periodic external input. 

Using the model equations ([l])-(l3|), we can uncover the 
bifurcation scenario which is responsible for the observed 
bistable behaviour. In Fig. [3] (a) and (b) the local ex- 
trema of the numerically calculated field intensities 
and I i?2 P are shown for increasing and decreasing K and 
fixed Aoj. The bifurcations evident in Fig. |4](a) and (b) 
can be understood from the global bifurcation diagram 
in the Aui vs. K plane shown in Fig. [4] (c), which was 
calculated using the numerical continuation tool AUTO- 
07p ^21!. As K increases from zero, the stable two mode 
limit cycle associated with Alu which exists for small K 
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FIG. 4: Local extrema of the field intensities (a) \E\^ and 

(b) \E2^ obtained from numerical integration of Eqs. 

((31) as a function of the injected field strength for increasing 
and decreasing injection. The frequency detuning is -14 GHz. 

(c) Bifurcation diagram in the Alj vs. K plane. Subcritical 
and supercritical Hopf bifurcations are indicated by dashed 
and solid blue lines and labeled via and , respectively. 

and denote Hopf bifurcations of single and two mode 
equilibria, respectively. The solid black line SN is the saddle- 
node bifurcation of the single mode locked state. 



evolves into a stable two mode equilibrium via a super- 
critical Hopf bifurcation at H'^ ■ With increasing K this 
two mode equilibrium becomes unstable via the subcriti- 
cal Hopf bifurcation at EI2 ■ On the other hand, however, 
the single mode locked state becomes stable via a saddle- 
node bifurcation at the fine SN (cf. Fig. [4](c)). Bista- 
bility therefore occurs if the SN bifurcation occurs at a 
smaller K value than the fl^ bifurcation, which is the 
case for sufficiently negative detuning (Aw < — lOGHz) 
as evident from Figd] (c). This explains the hysteresis 
loops found in Fig. [4] (a) and (b) . 

We are now in a position to connect our numerical 
findings in Fig. U] with the experimental results of Fig. 
[5] For sufficiently large Aw, the four eigenvalues of the 
Jacobian matrix of the two-mode equilibrium comprise 
two complex conjugate pairs with imaginary parts equal 
to the detuning frequency Aw and the relaxation oscilla- 
tion frequency z/ro ■ We identify the boundary K\ in Fig. 
[5] with the supercritical Hopf bifurcation Ei^ in Fig. U) 
At the real parts of the eigenvalues associated with 
Aw change sign and become negative with increasing K. 
Experimentally this is reflected in a signiflcant broaden- 
ing of the spectral feature at frequency Aw in the dy- 
namics of the injected mode. This is consistent with the 
expected behaviour in a noisy system and explains the 
observed double peak structure in the optical spectrum 
shown in the upper left inset of Fig. [TJ We identify the 
upper boundary of the bistable region K-^ in Fig. [2] with 
the subcritical Hopf bifurcation of Fig. U) At the 
real parts of the eigenvalues associated with v^o change 
sign. Within the bistable region, spontaneous emission 
noise is able to excite precursor oscillations at fpo which 



are visible in the power spectra of both modes in Fig. [2l 
Finally, the lower end of the bistable region K2 in Fig. [2] 
can be identified with the saddle-node bifurcation SN in 
Fig. m When approaching this bifurcation with decreas- 
ing no precursor oscillations are observed, as saddle- 
node bifurcations do not generically involve oscillatory 
behaviour. All of the above conclusions are confirmed by 
numerical simulations that include spontaneous emission 
noise. 

It is interesting to compare this bifurcation structure 
with others considered in the literature. In [l^ a bista- 
bility between a limit cycle and a locked state was re- 
ported in optically injected single mode lasers. This limit 
cycle becomes unstable via a subcritical torus bifurca- 
tion. In contrast, our bistability mechanism involves only 
Hopf and saddle-node bifurcations and occurs between 
two equilibria. The relative simplicity of our bifurcation 
scenario is therefore a result of the presence of the sec- 
ond lasing mode, and the increase of the dimensionality 
of the dynamical system by one. A pair of subcritical 
and supercritical Hopf bifurcations have also been asso- 
ciated with relaxation-oscillation and wave-mixing pul- 
sations in the context of in-plane semiconductor lasers 
with integrated optical feedback [22| where similar pre- 
cursor features close to the bifurcation points were also 
reported. 

The bifurcations of a two-mode system have also been 
studied in the context of optically injected VCSELs to 
explain the switching between single-mode states with 
different polarizations 



14| 



However, since the two po- 
larization modes in a VCSEL have a very small frequency 
spacing (c. 10 GHz), these polarization modes are intrin- 
sically bistable, and a six dimensional system of equa- 
tions is required for adequate modeling. As a result the 
bifurcation scenario differs significantly from that consid- 
ered here, which involves an equilibrium state with two 
simultaneously active and highly nondegenerate lasing 
modes. We note also that the case of two-mode injection 
locking in an optically injected VCSEL [2^ must also 
be contrasted with our case, as injection locking implies 
a frequency degeneracy between the orthogonally polar- 
ized modes and the injected field. Such a state, which 
can be regarded as a single mode state with elliptical po- 
larization, cannot arise in the context of our model as the 
phase of the uninjected mode is decoupled. 

In Fig. [5] we show the results of numerical modeling of 
the optical memory element. For this simulation we have 
used a 5 ns pulse duration and we also included spon- 
taneous emission noise. Agreement with the measured 
intensity time traces is very good. In particular, we have 
reproduced the large and decaying relaxation oscillations 
in the injected mode's unlocking dynamics. The fast lock- 
ing of the injected mode and the fast extinction of the 
uninjected mode is also reproduced as is the rather slower 
recovery of the uninjected mode once the injected mode 
unlocks. Numerical results suggest that the memory ele- 
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FIG. 5: Numerical intensity time traces of the injected field 
strength and of the two primary modes of the device. The 
detuning is -14 GHz. The average value of the injected field 
strength is 0.018 with a modulation of ± 0.008. The pulses 
are 5 ns duration with a 200 ps rise-time. 



ment can be switched between states with subnanosecond 
pulses. 



CONCLUSIONS 

We have demonstrated an all-optical memory element 
based on the injection locking of a two-color Fabry-Perot 
laser diode. Pulsed modulation of the injected field leads 
to switching of the uninjected mode with 35 dB intensity 
contrast. The underlying bistability in the system was 
found to be between the single mode injection locked 
state and a two-color equilibrium state, unlike in the 
single mode injected laser. Numerical modeling of the 
switching dynamics was in excellent agreement with ex- 
perimental results. The simplicity of our approach should 
lead to new opportunities for all-optical signal processing 
with in-plane semiconductor lasers. 
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